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ABSTRACT. The 1984 article by Martin, Odlyzko and Wolfram on Wolfram’s Rule 
90 heightened interest in the area of linear cellular automata defined on a vector 
space over a finite field [6]. In this article we find the closed form expression for 
the minimal polynomial of Rule 90, which is then used to completely characterize 
the dynamics of the map. In addition, we show that the standard basis vectors in 
Z3 lead to the maximal cycle for any linear finite dimensional cellular automata 
with periodic boundary conditions. Finally, we address questions posed by Tadaki 


on the connection between Rule 90 and Rule 150. 


1. INTRODUCTION 


A finite dimensional cellular automaton (CA) is a discrete time dynamical system 
defined on a finite dimensional vector space for which the next state is updated by a 
local deterministic rule. Each component in the vector space is considered as a cell 
typically taking on only two values, 0 and 1. If the vector space is defined over a finite 


field, all forward orbits of a CA converge to a periodic cycle in finite time. Therefore, 
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for finite dimensional CA, the main problem under consideration is determining the 
transient behavior and cycle lengths of the map. 

In the case of a CA defined over a finite dimensional string taking on values from 
{0,1}, an important distinction is whether the update map acts linearly when the 
state space is viewed as a vector space over Zp. In such a setting, one can exploit 
the algebraic structure of a vector space over a finite field and apply the many known 
results of iteration of linear maps. 

The concept of CA as a dynamical system on a vector space over a finite field was 
employed by Martin, Odlyzko and Wolfram in their 1984 paper [6] on Wolfram’s Rule 


90 given by: 


WrX = (fn + T2, 21+ 23,..-,%n-1 + 21) (1) 


where x = (z1, %2,...,£n) E€ Z} and the addition is modulo 2. Since 1984, many 
articles have appeared extending and generalizing the results in [6] to other linear rules 
[4,5,8,9]. The focus of the research on W,, and its extensions involved characterizing 
and predicting the cycle lengths of the linear CA by utilizing the finite field structure. 

E. Jen characterized limit cycle structure through orders of minimal polynomials 
in the cylindrical case. Stevens proved that all cycle length for any finite dimen- 
sional linear CA can be obtained as orders of minimal annihilating polynomials [9]. 
Moreover, the transient behavior can be obtained from the number of factors of A 


belonging to the minimal polynomial. 


In two articles [10,11], Tadaki uses a recursion equation to find the characteristic 
polynomial of Rule 150 and connect the cyclic behavior to Rule 90. However, as 
stated earlier, period lengths are not obtained through the characteristic polynomial, 
but through the minimal polynomial. To this end, this paper restates the proof that 
all period lengths are equal to orders of factors of the minimal polynomial of the 
map. In Section 3, we find a closed form expression for the minimal polynomial of 
Rule 90 using a dimensional and symmetry argument. Several observations on Rule 
90 are direct results of the proof of this theorem. Section 4 discusses Rule 150 and 
its connection to Rule 90. Data on all periods up to n = 40 for Rules 90 and 150 are 


provided. Some general problems for further work are suggested in Section 5. 


2. CHARACTERIZATION OF CYCLE LENGTHS AS ORDERS OF MINIMAL 


ANNIHILATING POLYNOMIALS 


Although the results in this section are posed over Z} due to the present appli- 
cations, all results hold for any generalized finite field with characteristic p. The 


following definitions are used extensively in our analysis. 


Definition 2.1. The minimal annihilating polynomial of a vector v € Z3 is the 


monic polynomial j1,(A) of least degree such that jz,(A)v = 0. 


The existence of such a polynomial is guaranteed by the Cayley-Hamilton theorem 


which states that the characteristic polynomial of A will annihilate A. 
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Definition 2.2. Let q(A) be a polynomial with coefficients in Z2. Suppose that 
q(0) # 0. Then the order of q(A), ord(q(A)), is defined to be the smallest natural 
number, c, such that q(A)|A°— 1. If q(0) = 0, then q(A) can be written as \*q(A), for 
some positive integer k, where the polynomial, g(A), has the property, g(0) 4 0. In 


this case, the order of g(A) is defined to be the order of q(A). 


The characterization result that identifies transient behavior and cycle lengths fol- 


lows: 


Theorem 2.1. Let v € Z}. Let (à) be the minimal annihilating polynomial of 
v. Assume that (A) = AFT (A) where k > 0 and (LAT is a monic polynomial with 
LUIT A 0. Then the kth iterate of v belongs to a periodic cycle with period length 


c= ord( mw). 


A variation of this characterization result first appeared in [5]. We restate the proof 
of this theorem that can be found in [9]. 
Proof: 
Let A’v be the first iterate that belongs to the periodic cycle. Denote the minimal 
length of the cycle by c. Then A°(Av) = Atv, so A?(A°—I)v = 0. Since the minimal 
annihilating polynomial divides any other annihilating polynomial of v we know that, 
Hol AAO = 1). 


This means that ANI and 7,(A)|\° — 1. Therefore, ord(i,(A)) < c. 
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Now we will show that ord(/i,(A)) = c. Assume on the contrary that ord(/i,(A)) = 
L< c. Then ji,(A)|\'—1, A*®(\!— 1) = py (A)q(A) for some q. Therefore A*(A! I)v = 
Ly(A)q(A)v = q(A)p,(A)(v) = 0, which means that A’A*v = A*v. Therefore, A*v 
is on a periodic cycle of length l < c. This is a contradiction to the minimality of c. 

Now we will show that k = j. Recall that Afv was first iterate on cycle and 
HLA) = A*fio(A). 

Because A*|)\ it follows that k < j. Assume on the contrary that k < j. 

Given that ,i,(A)|A° — 1, we have py(A)|A*(A° — 1). This implies ARLA" — 1) = 
[y(A)p(A), for some polynomial p. Therefore, A*(A°—1)v = p(A),(A)v = 0. Hence, 
ALAR) = A*v, which implies that A*v is on the periodic cycle. Recall that A’ (v) 


is first iterate to be on the cycle and so we arrive at a contradiction. a 


The minimal annihilating polynomial for any given v is a factor of the minimal 
polynomial. Therefore, all possible period lengths can be obtained from the minimal 
polynomial of A. Moreover, the maximal period length is equal to the order of 
the minimal polynomial because there exists a vector whose minimal annihilating 
polynomial is equal to the minimal polynomial [2]. 

It is important from a dynamical systems perspective to connect the iterates of a 
vector v under A to the minimal annihilating polynomial of v. If v € ZY is a vector 
contained in a cycle under A, then the algebraic period of v under A is the least 


positive integer value a such that Av is a linear combination of v, Av, A2v,...A%!v. 
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If a is the algebraic period of v then we may write 


Atv = bu + ba Av K... ba Anh, 


Then the minimal annihilating polynomial of A is given by 


N= bg-1A%* =... hah h 


Thus, a equals the degree of j1,(A). The concept of algebraic period was discussed in 
[3], although not named as such. 


We now use the characterization result to analyze the CA of interest. 


3. WOLFRAM’S RULE 90 


Wolfram’s Rule 90, defined by the map (1) was first studied in the finite field setting 
in [6]. The map first appears in [7] as a model to predict stunted tree growth in a 
forest. We now discuss W, using the language of minimal polynomials. 


The matrix representation of W, in the standard basis is, 


01 0 01 
10 1 0 0 

Awn = (2) 
0 101 
1 0 Br L 0 


An algorithm based on Theorem 2.1 was developed in [9] to generate all cycle lengths 
of a linear CA. Based on this algorithm applied to Awn cycle lengths up to n = 40 


were computed and appear in Table 4. 


3.1. The Characteristic Polynomial for Wolfram’s Rule 90. 


3.1.1. The Null Boundary Map. Wolfram’s Rule 90 can be thought of as having peri- 
odic boundary conditions or as an infinite sequence on a cylinder. Maps on cylinders 
have been extensively studied in the literature as in [1,4,5,6,8,9,13]. In an effort to 
generalize the minimal polynomial result to maps without periodic boundary condi- 


tions, Stevens et al. [8, 9] examined the Null Boundary Map, defined by: 


Na (x) = (£2, £1 + £2, 22 + L3 + . . . , En—2 + Ln, Ln-1) 


where x = (£1, £2, Tal € Z3. 
The purpose of looking at N, in this paper is to utilize relationships between the 
characteristic polynomials of Ay,,, and Aw, to obtain a closed form expression for the 


characteristic polynomial of Ay,,. The following theorem illustrates this connection. 


Theorem 3.1. Let pw alL) be the characteristic polynomial for the nx n matrix Aw 
and pyn—1(A) be the characteristic polynomial for the (n—1) x (n—1) matriz Ann-1. 


Then pwn(A) = APN,n-1 (A). 


Proof: 


To compute the characteristic polynomial of Aw,, we need to calculate the determi- 


nant, 


A 1 0 0 0 
Toe 1 0 0 

|Awn T AMI = d : (3) 
0 1 Ad 
1 0 Oro Ds 


The determinant expands as 


M NAM Ma + |Minl 


where M;; is the n — 1 x n— 1 minor obtained by eliminating row 7 and column j 
from the matrix Awn + Al. 
Direct inspection reveals AM 11 = Apn n-1(A). 


The transpose of Mj is the matrix 


X 1 0 0 1 

1A 1 0 0 
Mi= 

0 L X 1 

0 0 01 Â 


Since, interchanging rows does not change the value of the determinant, we have that 


1A 1 U 0 1 0 U 0 1 
1 A 1 U 0 
Min = Jo sa Eh oe WS 
0 0 O 1 A 0 1 à 1 
10 0 ... 0 1 U 0... O 1 A 
= IAM = |M.. 
Because the computations are over Z2, |Mi 2| + |Min| = 2|Mi2| = 0. E 


The characteristic polynomial for Ay „n was obtained in [8] by applying the method 


of generating functions to the recursion formula: 


PN n LAT = AP Nyn — PN,n-1- 


The closed form expression for the characteristic polynomial is given by, 


i= es a) 


[n/2] 


where all the binomial coefficients are to be interpreted modulo 2. 


Theorem 3.2. Letn = 2k. Then the characteristic polynomial for Awn is given by 


pwn(A) = Ces + (i+ OPLEES Gs ey (5) 
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Ifn =2k +1 then 


Proof: 
The closed form for py,,, is found by obtaining the characteristic polynomial through 
the relationship formulated in Theorem 3.1 and substituting n = 2k, 2k + 1 into (4). 


The substitution and repeated applications of Kummer’s theorem yields, 


Wee Cary eM Ge ha eee a ifn =2k+1 
Dun LA) — 
(A (Ges) =o (FJA T es w+... + eC AH) if n is even. 


Since the minimal polynomial is a factor of the characteristic polynomial, we will use 
the characteristic polynomial to arrive at the closed form expression for the minimal 


polynomial. 


3.2. The minimal polynomial for Wolfram’s Rule 90. The following theorem 


gives the minimal polynomial for Ay». 
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Theorem 3.3. Let mwa LA) be the minimal polynomial for the n x n matrix Awn. 


Let pw. be the characteristic polynomial for the n x n matriz Awn. Then 


Pwr) if n = 2k; 
Lwal(A) = (7) 


APN lA) — Ppnwr(A)) ifn = 2k +1. 
Proof: 
Case 1: n = 2k 
Step 1 Obtain the relationship between n = 2k and n=k 
As stated in Section 3 we know that pywn(A) = Apnn—i(A). Furthermore, it was 


proved in [8] that py,2.-1 = Apy.,_1- Substituting yields, 


Pwk(A) = APN e-l) = Pw (A): 


Step 2 Verify the structure of the Smith Normal Form of Awn 
Because Aw, — AJ has an n — 2 x n — 2 minor with determinant one (Figure 2), we 


know that the Smith Normal Form has the structure 


1 0 0 0 0 
0 1 0 0 0 
Sn = (8) 
1 0 0 
0 0 PALA) U 
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where rn(A)Sn(A) = Pw lA), Tr(A)|Sn(A) and Sn(A) = mw» LA, 
Step 3 Show deg(ra(à)) = deg(sn(à)) < k. 
We will illustrate the concept of our argument by discussing the problem on ZY first. 
Consider the basis vector e; = (1,0,0,0,0,0). Its components can be assigned to 
the vertices of the polygon having the corresponding number of vertices, as shown in 


Figure 1. The iterates of e} under Aw, are also depicted there. Because Awn is 


S) Ale J. ls) JD fe) 


FIGURE 1. Iterates of e; = (1,0, 0,0, 0,0). 


symmetric, all iterates are symmetric as is evident in Figure 1. Thus, in general, the 
maximum number of linearly independent iterates of e, is k and hence the algebraic 
period of e; is less than or equal to k. Thus, the degree of the minimal annihilating 
polynomial for e; is less than or equal to k. 

We claim that pe, LA) is the minimal polynomial of Awn. To see this first observe 


that pe, (A) = He, (À) = --- Me, (A). This is true because Aw, commutes with the shift 
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FIGURE 2. Shaded region is a n — 2 x n — 2 submatrix of Awn — AI 


with determinant one. 


map: 


Hei (Aw,n)ei = He, (Awn) Ste = S$ He, Lm ale =0. 


Moreover, if v € Z} then v = X; ce; where c; € Z3 and pe, (Awn)v = Xi] Cile, (Awn)ei = 
0. Therefore deg(uwn) < k. This proves that deg(s„) < k. The divisibility condition, 

Tn|Sn, implies the degree of r, is also less than or equal to k. 

Step 4 Show that r,(A) = sn(\), which proves that pwn(À) = a LA, 

Since r,(A)Sn(A) = pwn(A), we know that deg(r,) + deg(s,) = 2k. But since 
deg(rn) < deg(s,) < k, this implies that deg(r,) = deg(s,) = k. Since r,|s,, we 

have that rn(A) = sn(A). Hence, pwn(A) = Haa (A) = (pwx(A))*. This proves that 

HW (À) = pw,r(A). 

Case 2: n=2k+1 


The n odd case is proved in a similar fashion as the n even case with minor changes. 
13 


If n = 2k + 1 then from a result in [8], 
Pw,2k+1 = APN 2l A) = Alpu el A) + pN a)”. 


Let qn (A) = Pnl À) + HKL). 

Step 1 Show «LAT has only one factor of à. 

By the closed form expression, (4), py2x(0) = 1 and pyn 2k+1(0) = 0. Thus, qn(0) = 1. 
Therefore, \|pwn(A), but A? pwn(à). It follows that Aļsn (A) but A?)/s,(A), where 
Sn(A) is the (n—1,n-— 1) entry in the Smith Normal Form, (8). Write s,,(A) = AS, (A). 
Step 2 Show dog aw LAT <k +1. 

The proof as before uses symmetry arguments to show that deg(ue LAT < k +1. To 
see this we need to think of e; on a polygon as illustrated in the case n = 5 (Figure 3). 


Due to the symmetry, there can be only at most k +1 linearly independent iterates of 


l U (| U 
O U } ( i ( j 
U U U 1 1 1 ji 
S) Asla) Asla) Asle) 


FIGURE 3. Iterates of e; = (1,0, 0,0, 0,0). 


e1. Thus the dimension of span{e1, Aw.ne:,...} must be less than or equal to k + 1. 


Therefore, deg(ue, (A)) < k +1. As in Case 1, pe, is the minimal polynomial of Awn 
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and therefore deg (tw 2r41) < k +1. 

Step 3 Show pyn(A) = A(Sn(A))?. 

We know that deg(r,(A)) < deg(s,(A)) < k. We also know that deg(r,(A)) + 
deg(s,(A)) + 1 = 2k + 1. This implies deg(r,,(A)) = deg(S,(A)) = k. So we can 
conclude that rp(A) = §,(A) and hence, pw.n(A) = A(s2(A)) where Sv LA) = qn (A). 


A nice result of the above method of proof is that the standard basis vectors in Z3} 
lead to the maximum cycle for any n. In fact, the basis vectors in Z3 lead to the 
maximal cycle for any finite dimensional linear CA with periodic boundary conditions. 


To see this consider the general finite dimensional linear map with periodic boundary 


conditions: 
Ly (x) = (a1£1+a2£2+. . .+0n£1, Lot... +n -1En,..-,A201+...AnEn-1+G1%n) (9) 
where a; € Z fori = 1,...,n. The matrix representation of Ln in the standard basis 


vectors is the right circulant matrix, 


a, ag a3... An 

An A, Q2 a3... An—1 
ALn = 

a2 a3... An ay 


Multiplying Az n and Ag, where Ag, is the matrix representation of the left shift 


map verifies that Azn commutes with Ag,. Commutation with the shift map was 
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the only criteria required to show that e;,2 = 1,...,n led to the maximal cycle under 
Wp. Thus, the same argument verifies that the basis vectors lead to the maximal 


cycle under L and therefore we can state the following theorem. 


Theorem 3.4. Let L, be the general linear CA defined in (9). Then the standard 


basis vectors, e;, i =1,...,n lead to the maximal cycle under Ly. 


We now examine Rule 150 and its connection to Wp. 


4. WOLFRAM’S RULE 150 AND THE BANDED MAP 


Wolfram’s Rule 150, 


T(x) = (En +21, +%2,%1 +%2+%3,...,€n-1 + Xn + zı) (10) 


where x = (£1, £2,... 2n) € Z was studied by Tadaki in [10,11]. In [11], Tadaki 
arrived at a recursion formula for the characteristic polynomial of Arn. In actuality, 


Tadaki was working with the null boundary condition map 


Bp (xX) = (£1 + £2, £1 + £2 + T3, . . Daa ag + Pais Pal + Tn), 


although not formally stated. 

Tadaki’s goal was to predict cyclic behavior of Rule 150 by using the closed form 
expression for the characteristic equation and then to find the relationship between the 
cycle lengths of Rule 90 and Rule 150. As proved in Theorem 2.1, cyclic behavior is 


characterized through the minimal polynomial and thus a closed form for the minimal 
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polynomial of Arn is desired. Because Arp, = uw aY], the characteristic and minimal 
polynomials of Rule 150 can be obtained from pus LA) and uwn(à). In fact, we have 
that prn(A) = |Awntl—-Al| = |Awn-(A-1)L = pwn(Atl). The same substitution 
yields prn(A) = Lwn(A+1). This relationship leads to the following result connecting 


the behavior of both maps. 


Theorem 4.1. [fn = 2", the order of twrn(A) is 2471. Furthermore, the order of any 


minimal annihilating polynomial is 27 where j < k — 1. 


Proof: 


It is shown in [8] that 
Pn 2k-1(A) = aes 
Therefore, 
Pw2k = ae 


By the closed form expression for the minimal polynomial of Awn given by, (7), 


k-1 
Hw, = xe 


gk-1 


Thus ord(uwn(à)) = 1. From this formulation we have that prn(A) = (A+ 1)". 
The equality implies that ordi «LAT = 2-1. From the closed form expression, we 


also can see that (\ + 1)” is a factor of prn(à) for j < k — 1. Consequently, all 
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Vector Length | Cycle Lengths | Cycle Lengths || Vector Length | Cycle Lengths | Cycle Lengths 


Rule 150 Rule 90 Rule 150 Rule 90 


R HT R 
RE HE E T G 
SEO O a GE 
HH | e | ae | ea 
a ae HE E e a 
Cie [ase | ee Perm | s s 
Cen | o | am S a enea aeaa 
Sa e HE e ae 718 
Tass | a | e [eee H | 
Tare [ans [an [eee [om | a 
Caas ee | os | n T 
Cae | as | ee [aaraa 
Cae | E | ea [ae [iaaa 
Tons | om | ram | nes [var [en 
Tas | ea | a [eee | nma [en 
Ceea [na | aen [oem [o [neve 
ea [ane | e G G aea 


TABLE 1. Period lengths under iterations of the Wolfram’s Rule 150 


and 90. 


minimal annihilating polynomials for n = 2* must have order 27 where j < k — 1. 
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5. CONCLUSION 


In this paper, we have obtained closed form expressions for the minimal polynomials 
of Rule 90 and Rule 150. Since cycle lengths and the transient dynamics can be 
obtained from the minimal polynomials, such expressions are valuable. Through the 
computation of the minimal polynomial of Rule 90, we find that the standard basis 
vectors always lead to the maximal cycle for any finite dimensional linear CA with 
periodic boundary conditions. 

Using the relationship that fy»(A+1) = UT nlà) we obtained connections between 
the period lengths for both maps when n is a power of two. It would be interesting 
to investigate further the connection between the order of an arbitrary polynomial, 


q(A), and the order of q(A + 1). 
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